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Abstract 

In the study of electronic transport in nanostructures, electric current is commonly considered 
homogeneous along the sample. We use a method to calculate the magnetoresistance of magnetic 
nanostructures where current density may vary in space. The current distribution is calculated 
numerically by combining micromagnetic simulations with an associated resistor network and by 
solving the latter with a relaxation method. As an example, we consider a Permalloy disk exhibiting 
a vortex-like magnetization profile. We find that the current density is inhomogeneous along the 
disk, and that during the core magnetization reversal it is concentrated towards the center of the 
vortex and is repelled by the antivortex. We then consider the effects of the inhomogeneous current 
density on spin-torque transfer. The numerical value of the critical current density necessary to 
produce vortex core reversal is smaller than the one that do not take the inhomogeneity into 
account. 
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I. INTRODUCTION 



Electric transport in magnetic nanostructures is a useful tool both for probing and for ma- 
nipulating the magnetization. In the low current density regime, magnetoresistance curves 



i current density 
3J. Magnetore- 



are useful for probing the sample's magnetization state while, in the higl 
regime, magnetization patterns can be modified by a spin-transfer torque 
sistance measurements have the advantage of being relatively simple and fast, serving as an 
efficient magnetic reading mechanism 0,0]. 

Depending on thickness and diameter, small ferromagnetic disks exhibit stable topological 
defects called magnetic vortices These vortices can be manipulated by picosecond 

pulses of few (tens of) Oersted in-plane magnetic fields that switch their polarity js- 13], 
making them good candidates for elementary data storage units 9]. 

For their use as storage units, the most viable form of manipulation of the magnetization 
is through spin-torque transfer, with the injection of high density electrical currents 
The effect of these currents in the magnetization dynamics is described theoretically by 
the incorporation of adiabatic and non-adiabatic spin-torque terms in the Landau-Lifshitz- 
Gilbert (LLG) equation 14, [ijj]. These two terms are proportional to the injected current 
density and it is normally considered an homogeneous current distribution inside the disk. 
Although theoretical predictions using this approach agree qualitatively with experimental 
results, there is a lack of quantitative agreement between theoretical and experimental results 
regarding to the current densities necessary to modify the magnetic structures 17H19]. 

In this paper we investigate the effect of non-uniform current distributions on electronic 
transport and spin-torque transfer in ferromagnetic systems exhibiting vortices. We calculate 
numerically the magnetoresistance (MR) and local current distribution of a ferromagnetic 
disk by separating the timescales for magnetic ordering and electronic transport. We consider 
an effective anisotropic magnetoresistance (AMR) that depends on the local magnetization. 
We discretize the disk in cells and solve the Landau-Lifshitz-Gilbert (LLG) equation 20 ] 
numerically with fourth-order Runge-Kutta 21], obtaining the magnetization profile of the 
disk. This pattern is used to calculate the magnetoresistance of each cell as a fixed current 
I is applied at two symmetrically distributed electrical contacts, resulting in a voltage drop 
and an inhomogeneous current distribution along the disk. 

This method couples the electric and magnetic properties of the metallic nanomagnets 
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and can be used to analyze the effect of inhomogeneous current distributions in different con- 
texts. First, we discuss the limit of low current density where transport measurements can 
be used to probe the magnetic structure. We compare the magnetic structure with magne- 
toresistance curves and show how magnetoresistance measurements could be interpreted to 
obtain information on the magnetization profile and its dynamics during the vortex core mag- 
netization reversal. Moreover, we discuss the consequences of a non-homogeneous current 
distribution on spin-torque transfer and find that the critical current density that produces 
vortex core reversal is reduced by one order of magnitude whenever such non-inhomogeneity 
is taken into account. This result can be seen as a new route to understand why experimen- 
tal values of the critical current densities are usually lower than the ones obtained in LLG 



calculations 



17H19|. 



This article is organized as follows: In section [TT] we discuss the model and method 
for the calculation of the magnetoresistance and current distribution. In section 1111} we 
exemplify the calculations by considering the magnetoresistance and current distributions 
of a Permalloy disk exhibiting a magnetic vortex. In section ITVl we study the consequences of 
a non- homogeneous current distribution on spin-torque transfer. In section |V] we summarize 
the main results. 



II. MAGNETORESISTANCE AND CURRENT DISTRIBUTION CALCULA- 
TIONS 

Let us consider a 36nm-thick Permalloy disk with a diameter of 300 nm discretized into 
a grid of 4 x 4 x 4nm 3 cells. The dynamics of the magnetization vector associated with 
each cell is given by the Landau-Lifshitz-Gilbert equation, which we numerically integrate 
with fourth-order Runge-Kutta and discretization step h = 10~ 4 {2]]]. The parameters 
associated with the LLG equation are the saturation magnetization M s = 8.6 x 10 5 A/m, 
exchange coupling A = 1.3 x 10 _11 J/m and Gilbert damping constant a = 0.05 [13j. 

By varying the external in-plane magnetic field H from negative to positive saturation 
we obtain a hysteresis curve, as depicted in Fig. [TJ which is consistent with experimental 
observations [7]. As shown in Fig. HJa), in static equilibrium and in the absence of magnetic 
fields, a vortex structure with a core magnetized perpendicular to the disk plane is formed in 
the center of the disk. If a small in-plane magnetic field H is applied, the core is displaced 



3 




4 



time scale of electronic scattering processes. 




R.=R X 



FIG. 2: Original cells used in the LLG simulation with the associated resistance network. 



Guided by recent experiments |23| we allow a constant current / to flow along the disk 
by attaching symmetrically placed electrodes on it (see Fig. [3]). The voltage drop along 
the resistors and the associated current map of the disk are obtained by solving Kirchoff 's 
equation iteratively at each node of the grid with a relaxation method 24|, |25] : 



V; 



(n+1) 





(1) 



where Rij is R^(R^) if i and j are horizontal (vertical) neighbors and 6, is the boundary 
current, assumed to be at the leftmost (rightmost) cells and zero otherwise (see Fig. |2]). 

V™ is the voltage at site % after n iterations and the sums run over the nearest-neighbors 
(j) of node i. Starting with a random initial condition {V^ -*} at each site we iterate EqffJ 
until each becomes stationary (within 9 decimal digits precision). After convergence 
we calculate the equivalent resistance, the ratio R eq = AV/I between the voltage drop AV 
between the electrodes, given by 



AV = E y i- E v * 

i\\bi=I j\\bj=—I 

and the current I entering the disk. 



(2) 
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III. MAGNETO-STRUCTURE AND MAGNETORESISTANCE 



A. Hysteresis and magnetoresistance 

In order to obtain the magnetoresistance curves, the calculation discussed in the previous 
section is performed at different fields. Magnetoresistance and current distribution for the 
same points of the hysteresis loop in Fig. [1] are depicted in Fig. [3j Fig. |3]^a) displays 



magnetoresistance curves for both homogeneous (without using the resistance network 26]) 
and non-homogeneous current distributions. The vortex expulsion and its formation at a 
different critical field are clearly identified and, with p" = 155f2-nm and p 1 - = 150f2-nm, we 
obtain a MR of 1.2% for the non-homogeneous distribution, which is a typical value found 

1,0- 



m experiments 



23 



One also observes that magnetoresistance curves for uniform and non-uniform current 

y, the latter being more comparable to experimental results 



distributions differ significant 
with same contact geometry 23J. As expected, a homogeneous current overestimates the 
magnetoresistance, since the current will flow thru regions of high resistance, whereas if the 
current found by solving Laplace's equation on the associated resistor network, one finds 
a preferential path (higher current density) on regions of low resistance. This difference is 
more pronounced in the presence of a vortex, since the magnetization of the disk is highly 
non- homogeneous on such configuration. 

In the light of the discussion above, one sees in Fig. [3] (b)-(c) that the current is not 
homogeneously distributed inside the disk, being stronger towards the center of the vortex 
core. In the center of the disk the magnetization either points in z direction, perpendicular 
to the direction of current flow, or loops about the vortex core. In both cases, the current 
has a path where its direction is always perpendicular to the magnetization, reducing the 
local magnetoresistance. Above the saturation field, the magnetization is uniform and at 
the disk center the same happens to the current. The red (blue) region has a current density 
1 — 2% larger (smaller) than the current / at saturation. This effect might be enhanced if 
other sources of magnetoresistance are considered, as giant magnetoresistance for example. 
Similar approaches, using different sources of magnetoresistance and geometries had been 



used to calculate the magnetoresistance in nanomagnets [27H32||. 
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FIG. 3: (a) Magnetoresistance for the magnetic configurations obtained in Fig[T]for uniform (circle) 
and non-uniform (square) current distribution. Bottom: electric current map for (b) zero field and 
(c) H = 75 mT. The red (blue) color corresponds a current density which is about 1 — 2% larger 
(smaller) than the uniform current at the saturation field. 

B. Dynamics 



Next, we study the dynamics of the vortex core magnetization reversal by the application 
of short in-plane magnetic fields. Under a pulsed in-plane magnetic field or spin polarized 
current excitation, the vortex with a given polarity (V + ) dislocates from the center of the 
disk with nucleation of a vortex (V~)-antivortex (AV~) pair with opposite polarity after 
the vortex attains a critical velocity of rotation about the disk center [3|, |33[ . The original 
V + then annihilates with the AV~, and a vortex with reversed core magnetization (V~) {9] 
remains. If a low-density electronic current is made to flow through the sample (without 
disturbing the magnetization dynamics), we observe changes in the magnetoresistance, as 
the vortices nucleate and annihilate. We depict in Fig. H^a) the dynamics of the magnetore- 
sistance as a pulsed, in-plane magnetic field is applied in the x direction at t = 20 ps for 
different pulse intensities. The pulses have their shape sketched in grey in Fig. H^a) with 
full width at half maximum of t=250 ps. Depending on the pulse intensity, the vortex core 
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magnetization does not reverse at all (uqH < 43 mT), reverses once (54 mT< fioH < 64 
mT) or multiple times (noH > 64 mT) 12|, |l3j. 




FIG. 4: (a) Evolution of magnetoresistance after the application of pulsed, in-plane magnetic 
fields (the shape is shown in grey) with different intensities. Snapshots of magnetization (arrows) 
and current distribution (color map) for pulse fields without (b) and with (c)-(e) vortex core 
magnetization reversal. The white cross shows the position of the disk center. Both current and 
magnetic field are applied in the x direction 

During the application of a field pulse in the x direction, i.e., parallel to the current 
flow, the vortex core is pushed to the y direction breaking the rotation symmetry of the 
disk's magnetization, increasing both the total m x component and the disk's equivalent 
resistance (see Fig. H](a)). At t = 340 ps the field is practically zero and, from the decay of 
magnetoresistance to its equilibrium (initial) value, one can infer whether there was reversal 
of the vortex core polarization or not: for pulses that induce reversal, the value of the 
magnetoresistance just after the pulse is always larger than its initial value. If there is no 
reversal the magnetoresistance attains a minimum value that is lower than its initial value, 
i.e., before the application of the pulse, and oscillates about it. 



S 



In Figs. H)^b)-(e) we depict snapshots of current (color map) and magnetization (arrows) 
distributions at time steps marked with black dots in Fig. Ufa), in situations with or without 
vortex core magnetization reversal. Whenever the pulse decreases its intensity, the total m x 
component and the equivalent resistance of the disk follow the same pattern (although with 
some time delay), because the vortex core tends to return to the disk center, where m x = 0. 
Panel (b) shows the current distribution and magnetization at a moment corresponding to 
the minimum of the resistance curve, for a field intensity hqH = 27mT, for which there is no 
vortex core reversal. There is a large region with m y magnetization (and small m x ) in the 
center of the disk. This region, together with the vortex core, creates a low resistance path 
for the electronic current, decreasing the equivalent resistance towards a value below the 
equilibrium resistance. Panels (c), (d) and (e) show magnetization and current distributions 
at different moments of the vortex core magnetization reversal for a situation where there is 
a single reversal {[IqH = 64mT). In panel (c) we depict the current distribution at the exact 
moment of nucleation of the V~AV~ pair, the initial stage of vortex core magnetization 
reversal. Panel (d) shows the spin waves emitted just after the V + AV~ annihilation, a 
process that occurs with energy dissipation. Such energy loss drives the vortex core to the 
disk center along with some small oscillations, mainly due to reflections of spin waves at the 
edges of the disk. It turns out that the resistance follow equivalent behavior: it decreases 
towards the initial resistance value and remains always above it. Panel (e) shows the current 
distribution after the field pulse has vanished. As can be seen, the time dependent resistance 
curves can give us an indication of the vortex reversal process. 

Let us discuss in further details the interplay between magnetization pattern and current 
distribution. In Fig. [5] (a) we show a snapshot of the current distribution during the vortex 
core magnetization reversal process, with the V + and the V~-AV _ pair with negative polar- 
ity. As shown in Figs. [31 HI and [5](b) the current is pushed to the vortex core, where m x = 
and, consequently, the local resistance is minimum. With the nucleation of the AV~ vortex 
(Fig. [5] (a)), m x gets larger than zero around it, with m y — > 0. As current flows in the x 
direction, it is repelled from the antivortex core. 

In the latter analysis we considered a particular orientation of the AV. However, as can be 
seen in Figs. [5]^c) and (d), depending on their orientation, antivortices can either attract (in 
the first case) or repel currents (in the latter case). Vortices are rotation invariant, and alway 
attract current towards their cores. It is important to point out that this difference in current 
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distributions might have important consequences in the high-density current spin-torque 
transfer acting on either a vortex or an antivortex. For instance, although the inversion 
process through spin-torque for an AV is equivalent to the one for a V, we should expect 
different current densities in each one, since currents can only penetrate the AV core at a 
particular orientation. 




FIG. 5: (a)-(d) Snapshots of magnetization (arrows) and current distribution (color map) of 
(a) the vortex core magnetization reversal process at t = 214 ps, (b) a vortex and its associated 
current distribution, and (c)-(d) antivortices rotated by 45 degrees with respect to each other and 
their associated current distribution. Depending on the relative orientation of the antivortex it can 
either focus (c) or repel (d) currents away from the center of the core. 

IV. SPIN-TORQUE TRANSFER 

In this section, we discuss the consequences of inhomegenous currents in the spin-torque 
transfer. In order to check how the current distribution is incorporated in the spin-torque 
terms of the modified LLG equation, we need to review a few steps of their derivations. 
It is important to note that in our approach, the only source of non-homogeneous current 
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distribution is the anisotropic magnetoresistance, as discussed in section [TTJ All other effects 
are neglected. 

The itinerant electrons spin operator satisfies the continuity equation 

1{ S ) + V.(J) = - 1 -(([ S ,H))) (3) 

where J is the spin current operator. The Hamiltonian H is the s-d Hamiltonian {H s d = 
— J ex s ■ S), where s and S/S = — M/M s are the spins of itinerant and localized electrons, 
and J ex is the exchange coupling strength between them. We define spin current density 



J = (J) = -(gfi B P/eM 8 )j e (T) <g> M) 
spin density is given by m = (s) 



to calculate the spin-torque 
density. We obtain 
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, where j e (r) is the current density, and the electron 



JM. We use the same approximations, previously used 
, I5I ]. with the new ingredient of non-homogeneous current 



^m=^[M(V-j e (r)) + (j e (r)-V)M] 

dt eM (4) 

m x M, 

M s 

where, M is the matrix magnetization, g is the Lande factor splitting, /is is the Bohr 
magneton, P is the spin current polarization of the ferromagnet, e is the electron charge. 
From the continuity equation for charges, the term containing V • j e (r) is always zero, even 
if the current density is not constant. As discussed previously, the same divergent is used 
to determine the current distribution in section [III This expression is exactly the same 
expression obtained previously but with j e (r) in the second term of the right side of the 
equation varying with r. This current distribution is introduced at the modified LLG that 
considers spin-torque transfer. Therefore, we obtain a spin-torque transfer where the current 
distribution is not uniform. 

To consider spin-torque transfer effects we include adiabatic and non-adiabatic spin torque 
terms in the LLG equation, 

d TT d 

— m = -7 m x H e// + am x — m 

dt dt (5) 

-(u ■ V)m + /3m x [(u ■ V)m], 

where, m=M/M s is the normalized local magnetization, a is a phenomenological damping 
constant, 70 is the gyroscopic ratio, H e ff is the effective field, which is composed of the 
applied external field, the demagnetization field, the anisotropy field and the exchange field. 
The first term describes the precession of the normalized local magnetization about the 
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efective field. The second term describes the relaxation of the normalized local magnetization 
and (3 is a dimensionless parameter that describes the strength of the non-adiabatic term, 



which we consider to be 0.5 
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161 ] . The velocity u(r) = (gPfi B /2eM s )j e (r) is a vector 



pointing parallel to the direction of electron flow and j e (r) is calculated using the procedure 
discussed in section HT1 




FIG. 6: (a): Core polarity as a function of current density for homogeneous (squares) and non- 
homogeneous current distributions (with different MR), (b) and (c): magnetization profiles during 
the inversion process at the critical current for both non-homogeneous (b) and homogeneous (c) 
current distributions. The color map represent the out-of-plane magnetization, m z , and the arrows 
represent the in-plane component. 

To explain the importance of our assumption about the current distribution let us analyze 
the critical current density jg, the minimal current density needed to produce a vortex core 
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reversal. For this purpose, we simulated the magnetization dynamics of a system subjected 
to a DC current with the modified LLG equation (equation [5]). In Fig. [6ta) one sees 
the vortex core polarity as a function of current density j e . The different curves represent 
situations of homogeneous current (squares) and of three different values of AMR where 
the magnetoresistance ranges from 2% to 10%. Such AMR, as discussed in the previous 
sections, shape the degree of current inhomogeneity throughout the disk. One can see that 
the critical current density j c e in our model is 3 — 10% smaller than the one obtained for 
uniform currents. These results suggest a new route, together with the non-adiabatic term, 
to explain the discrepancy between experimental results and theoretical calculations of the 
critical current density jg. 

Our analysis might also have important technological implications, since we observe a 
(almost linear) correlation between the current density necessary to produce a core inversion 
and the anisotropic magnetoresistance of the material. Thus, by increasing the AMR of the 
sample, one can decrease the critical current jg, which is strongly desirable in memory devices 
for the sake of low energy comsumption and minimal heat waste. 

Panels (b) and (c) of Fig. [6] show the magnetic configurations for the moment just 
before the V + -AV~ annihilation, at the critical current density, for the model with non- 
homogeneous and homogeneous current distributions, respectively. In the case of inhomege- 
neous currents, the fact that a V(AV) attracts (repels) the current affects the velocity and 
separation distance of the V-AV pair during the vortex core reversal. As antivortices repel 
currents, the current density at their core is smaller, making them slower than vortices. 
As a result, after the nucleation of the V~-AV - pair their separation occurs faster than in 
the case where the current density in the center of a V or an AV is the same, as usually 
considered in micromagnetic simulations. 

V. CONCLUSIONS 

We performed a realistic calculation of magnetoresistance effects in magnetic nanostruc- 
tures subject that takes into account inhomogeneous current densities. For that purpose, 
we adapted a numerical relaxation scheme for the Laplace equation to the solution of the 
LLG equation for the magnetization profile along a Permalloy disk. Our results suggest 
that resistance measurements might be useful to probe the dynamics of vortex core mag- 
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netization reversal, induced by short in-plane magnetic pulses. Moreover, we note that the 
difference between current distributions close to vortices and anti-vortices have significant 
consequences for the spin-torque transfer effect. The inhomogeneous current distribution 
inside the magnet reduces substantially the critical current density necessary to produce a 
vortex core reversal. We conclude that materials with large anisotropic magnetoresistance 
need lower current densities to modify their magnetic structure, a much deasireable feature 
for most modern memory devices. 
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